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Abstract 






From 2005 to 2009, 25 observations of Cyg X-1 were performed with Suzaku, achiev- 
ing a total exposure of 446 ks. In all observations, the source was found in the low/hard 
(^ ■ State, while the 1.5-12.0 keV count rate of the All-Sky Monitor onboard RXTE varied by 

^ ; a factor of ~ 3. In each observation, the 10-60 keV HXD-PIN spectrum and the 60-400 

cn I keV HXD-GSO spectrum were fitted successfully by a thermal Comptonization model plus 

Pr ! reflection by a thick neutral material. As the soft X-ray intensity increased, the Compton 

y-parameter was found to decrease from 1.0 to 0.6, while the solid angle of reflection to 
increase by ~ 30%. Also conducted was timing analysis over a frequency range of 10~^- 
, . 10 Hz. As the source became brighter in soft X-rays, the characteristic frequency of hard 

^ • X-ray variation increased from 0.03 to 0.3 Hz, while the fractional hard X-ray variation 

integrated over 10~^-10~^ Hz decreased by a factor of ~ 5. The signals in the 60-200 
keV band were generally found to vary on shorter time scales than those in the 10-60 keV 
band. These spectral and timing results can be consistently interpreted by presuming that 
increases in the mass accretion rate cause the Comptonizing hot corona to shrink, while the 
optically-thick disk to intrude deeper therein. 

Key words: black hole physics — accretion — stars: individual (Cygnus X-1) — X- 
ray:binaries 



1. Introduction 

Black hole binaries (BHBs) are known to reside mainly in either of two spectral states, the 
low/hard state (LHS) and the high/soft state (HSS), which appear when mass accretion rate is lower 
and higher, respectively (McClintock & Remillard 2006). In the HSS, the accreting matter is consid- 
ered to release its energy predominantly in soft X-rays as a multi-temperature disk emission (Mitsuda 
et al. 1984, Makishima et al, 1986) from an optically thick and geometrically thin accretion disk, 
called standard disk (Shakura & Sunyaev 1973). In the LHS, in contrast, the spectrum becomes much 
harder, and is characterized by a flat power-law with a photon index of F ~ 1.5 and a high energy 
cutoff at ~100 keV . This has been interpreted as inverse Comptonization of some soft seed photons 
by hot Maxwellian electrons (e.g., Shapiro et al. 1976; Sunyaev & Triimper 1979) in a geomatri- 
cally thick and optically thin accretion flow, or "corona", in which advection may be operating (e.g., 
Ichimaru 1977; Narayan & Yi 1995). 

A scenario of the LHS, called "disk-corona" view (e.g., Liang et al. 1977), assumes that the 
system still harbors a standard disk, which provides the seed photons for Comptonization. Although 
the "disk-corona" model became a standard interpretation of the LHS of BHBs, the overall geometry 
of the disk and the corona remained poorly understood, including in particular, whether the geomet- 
rically thin disk extends to the last stable orbit (e.g.. Miller et al. 2006a,b; Rykoff et al. 2007) or is 
truncated at some larger radii (e.g., Ebisawa et al. 1996; Poutanen et al. 1997; Gierlinski et al. 2008; 
Nowak et al. 201 1). Equally unsettled is the interpretation of the rapid variation (e.g., Miyamoto et 
al. 1989) in terms of the disk-corona view. It was not clear, either, how these components evolve 
as the system approaches the HSS through the LHS. These limitations have been caused mainly by 
observational difficulties, to obtain high-quality spectra on short time scales, and over an energy range 
wide enough to determine simultaneously and precisely both the disk emission around ~1 keV and 
hard continuum shapes including the high energy cutoff energy of ~ 100 keV. 

In 2005 October, a broad-band Suzaku observation of Cyg X-1 in the LHS was conducted for 
an exposure of 17 ks (Makishima et al. 2008; hereafter Paper I). The obtained 0.7-300 keV spectrum 
has been reproduced successfully by assuming that the disk is truncated at ^ 15i?G, and it provides 
the corona with the seed photons. Here, Rq = GMbh/c^ is the gravitational radius, with G the 
gravitational constant, Mbh the black hole mass, and c the speed of light. The data also required a 
reflection hump with a solid angle of D./2tt ~ 0.4 and a mildly broadened iron line with a breadth 
of 0" ~ 1 keV; this is consistent with the disk truncation at ~ 15 Rq. The corona was inferred to be 
highly inhomogenious, expressed by a two-zone approximation invoking optical depths of ~ 0.4 and 
~ L5, and the same electron temperature of ~ 100 keV. In addition, variations on a time scale of 
1 s were studied in terms of Comptonized emission using intensity-sorted spectral analysis. It was 
then indicated that the X-ray intensity flares up, on time scales of ~ 1 s, when a larger fraction of 
the disk is covered by the corona, while the disk itself stays almost constant. In other words, the 
rapid X-ray flickering of Cyg X- 1 in the LHS is primarily generated by fluctuations of the number 



of disk photons intercepted by the corona. This view is consistent with the one, based on the almost 
energy-independent fractional rms variation spectrum, that the variation in the LHS is generated by 
fluctuations in the seed photon input (Gierlihski & Zdziarski 2005; Gierlihski et al. 2010). Actually, 
figure 8(b) of Paper I (high-phase to low-phase spectral ratio), which is almost equivalent to such 
rms spectra, is rather flat, except in ^ 1 keV where the ratio decreases due to the dominance of the 
disk emission (Yamada 2011). Thus, the results of Paper I reinforced the disk-corona picture, and 
quantified its details. 

Given this snap-shot result, our next task is to clarify long-term evolution of the disk-corona 
configuration, which is a key to understanding the transition from the LHS to the HSS. For this 
purpose, we analyzed 25 Suzaku data sets of Cyg X-1 acquired over 2005-2009. Among the obser- 
vations, the spectral and timing properties in hard X-rays (10-200 keV) showed significant changes, 
even though the source remained in the LHS throughout. In the present paper, we limit our analysis to 
the data obtained with the Hard X-ray Detector (HXD; Takahashi et al. 2007; Kokubun et al. 2007). 
Results with the X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) will be reported elsewhere, 
after the data are fully corrected for pile up effects. 

2. Observation and Data Reduction 

2.1. Suzaku observations 

So far, Cyg X-1 has been observed with Suzaku on 25 occasions from 2005 October 5 to 2009 
December 17, with a total exposure of 446 ks. The exposure of individual observations ranges from 3 
ks to 45 ks, with an average of 18 ks. Intervals between a pair of consecutive observations range from 
~2 days to ~1 year. A log of these observations is given in table 1. The optical axes of the XIS and 
the HXD differ slightly by 3.5', and the first observation was performed at the XIS nominal position, 
while the others at the HXD one. 

The data to be analyzed in the present paper were obtained from the HXD, which consists of 
Si PIN photo-diodes (HXD-PIN) covering 10-70 keV and GSO scintiUation counters (HXD-GSO) 
covering 50-600 keV. Sometimes the exposure is shorter than was aimed, due to increased electronics 
noise in HXD-PIN which reduced the live time. The exposures after correcting for dead times are 
shown in table 1 . Among them. Observation 1 is the one which was conducted in the Performance 
Verification phase by the Suzaku team, and already reported in Paper I. The utilized HEADAS version 
is 6.9, unless otherwise stated. 

Figure 1 is a hardness-intensity diagram of Cyg X-1, constructed from long-term 5.0-12.0 
keV and 1.5-3.0 keV data of the All Sky Monitor (ASM) onboard Rossi X-ray Timing Explorer. As 
represented by black dots in figure 1, the present Suzaku observations were all conducted when Cyg 
X-1 was in the typical LHS. Hereafter, we designate the 1.5-12.0 keV ASM intensity as Casm and 
use it as a measure of the soft X-ray intensity of Cyg X-1, because the XIS data still need extensive 
corrections for the pile up effects. This Casm is also considered as a rough measure of the mass 



accretion rate. 

2.2. Data Reduction and Background Subtraction 

We processed the 25 data sets in the same manner. Specifically, the PIN and GSO events were 
selected by criteria of elevation angle > 5°, cutoff rigidity > 6 GV, and 500 s after and 180 s before 
the South Atlantic Anomaly. The unscreened GSO data were reprocessed using hxdpi and hxdgrade 
to take into account the recent updates in the GSO calibration (Yamada et al. 201 1, in press), while 
the PIN data did not need reprocessing. 

Since the HXD has neither imaging capability nor an offset detector, the background contribu- 
tion to the on-source data must be estimated by simulations. The HXD background is dominated by 
the non-X-ray background, hereafter NXB. That of PIN was simulated by "PINUDLCUNIT" model 
(Fukazawa et al. 2009). Using the PIN upper-discriminator hit rate as a major NXB indicator, this 
model can reproduce the 10-70 keV NXB to an accuracy of 3% or better for an exposure longer than 
10 ks. The GSO NXB was simulated by "LCFITDT" model (Fukazawa et al. 2009), which tries to 
reproduce the GSO background in 32 energy bands by analyzing long-term (~1 month) light curves 
therein. This method achieves an accuracy of ~ 1% in background reproduction when the exposure 
is longer than 10 ks. The cosmic X-ray background (Boldt et al. 1987) can be ignored in our analysis, 
since Cyg X-1 is so bright an object. 

2.3. Light Curves of the HXD Data 

We have selected three typical observations, 3,10, and 15 in table 1, and show their two-band 
HXD light curves in figure 2 in the increasing order of the PIN count rate (which happens to be also 
the increasing order of the observation number). After subtracting the background and applying dead 
time corrections, the HXD-PIN (10-60 keV) and the HXD-GSO (60-200 keV) count rates in these 
observations were typically 20-50 cts s ^^ and 20-30 cts s ~^, respectively. The light curves reveal 
intensity variations by 20-30% on time scales of ^ 200 s, although those on shorter time scales are 
difficult to resolve in these plots. Hereafter, we use these three observations as the representatives of 
the data acquired at different source intensities. 

Background-subtracted and dead-time-corrected spectra of these three observations are shown 
in figure 3, in the response folded form. In the PIN band, the background subtracted signals in these 
three observations exceed the NXB by a factor of ~ 5, even at 60 keV. In the GSO band, the source 
signal is detected up to 300-400 keV, since it therein exceeds 3% of the background while the typical 
background uncertainty is ~ 1% (Fukazawa et al. 2009). 

3. Spectral Analysis 

3.1. The HXD Hardness Ratio 

As a first-cut characterization of the HXD spectra taken in the 25 observations, we calculated 
ratios of the 60-200 keV average count rate detected by HXD-GSO in each observation over that 









Table 1. 


An observation log of Cy] 


gX-1. 






No. 


Date 


MJD 


ObsID 


Exposure 
(ks) 


D.V.* 


Pointing 


epoch^ 


Casm 
(cts s"^) 


1 


2005/10/05 


53648 


100036010 


18.1 


2.1 


XIS 


1 


28.6 


2 


2006/10/30 


54038 


401059010 


27.7 


2.0 


HXD 


3 


21.0 


3 


2007/04/30 


54220 


402072010 


40.2 


2.1 


HXD 


3 


14.9 


4 


2007/05/17 


54237 


402072020 


32.6 


2.0 


HXD 


3 


13.5 


5 


2008/04/18 


54574 


403065010 


29.0 


2.2 


HXD 


4 


15.3 


6 


2009/04/03 


54924 


404075010 


15.5 


2.3 


HXD 


5 


21.7 


7 


2009/04/08 


54929 


404075020 


13.1 


2.3 


HXD 


5 


21.8 


8 


2009/04/14 


54935 


404075030 


13.0 


2.3 


HXD 


5 


17.6 


9 


2009/04/23 


54944 


404075040 


16.3 


2.3 


HXD 


5 


19.7 


10 


2009/04/28 


54949 


404075050 


12.4 


2.3 


HXD 


5 


20.8 


11 


2009/05/06 


54957 


404075060 


15.3 


2.3 


HXD 


5 


22.4 


12 


2009/05/20 


54971 


404075080 


17.5 


2.4 


HXD 


5 


23.6 


13 


2009/05/25 


54976 


404075090 


16.1 


2.4 


HXD 


5 


31.8 


14 


2009/05/29 


54980 


404075100 


25.7 


2.4 


HXD 


5 


33.7 


15 


2009/06/02 


54984 


404075110 


15.0 


2.4 


HXD 


5 


45.2 


16 


2009/06/04 


54986 


404075120 


6.8 


2.4 


HXD 


5 


41.8 


17 


2009/10/21 


55125 


404075130 


15.2 


2.4 


HXD 


6 


20.0 


18 


2009/10/26 


55130 


404075140 


20.0 


2.4 


HXD 


6 


16.1 


19 


2009/11/03 


55138 


404075150 


14.1 


2.4 


HXD 


6 


16.6 


20 


2009/11/10 


55145 


404075160 


18.3 


2.4 


HXD 


6 


15.1 


21 


2009/11/17 


55152 


404075170 


18.7 


2.4 


HXD 


6 


20.2 


22 


2009/11/24 


55159 


404075180 


11.4 


2.4 


HXD 


6 


21.8 


23 


2009/12/01 


55166 


404075190 


13.2 


2.4 


HXD 


6 


25.1 


24 


2009/12/08 


55173 


404075200 


17.4 


2.4 


HXD 


6 


22.8 


25 


2009/12/17 


55182 


404075070 


0.3 


2.4 


HXD 


6 


40.1 



D.V. stands for data version of the HXD data which differ in caUbration and data process. 
' Epoch is a period in which the same PIN response is applied. 



in the 10-60 keV by HXD-PIN, with the NXB subtracted from both. This quantity, called HXD 
hardness ratio, can be used as an indicator of spectral hardness in the hard X-ray band. 

As figure 4 shows, the obtained hardness ratio is clearly anti-correlated with Casm; this is 
mainly due to negative Casm dependence of the PIN count rate, while the GSO count rate is rather 
constant. This result extends the intensity versus hardness anti-correlation seen in the ASM range 
(figure 1) to the HXD band. 

3.2. Spectral Fitting 

For detailed spectral analysis, we prepared response files for HXD-PIN and HXD-GSO us- 
ing xspecl2. As listed in table 1, the HXD-PIN responses were selected according to "epochs" 
of the observation, which are updated when lower-discriminator levels and/or high-voltage val- 
ues of HXD-PIN are altered (Takahashi et al. 2007; Kokubun et al. 2007). The source po- 
sitions in the field of view, explained in subsection 2.1, also affect the responses. We applied 
ae_hxd_pinXXnomeY_20100107.rsp to the PIN spectra and aeJixd_gsoXXnom_20110107.rsp and 
ae_hxd_gsoXXnom_crab_20100526.arf to those of HXD-GSO. Here, XX and Y specify the source 
positions (XIS nominal as "xi" and HXD nominal as "hx") and the epoch ("1" to "6"), respectively. 
The arf file for HXD-GSO is a correction file to adjust the normalization between HXD-PIN and 
HXD-GSO, and to trim high energy responses (Takahashi et al. 2007). The PIN data below 14 keV 
were discarded since electrical noises can contaminate the data up to this energy. The GSO data below 
70 keV and above 400 keV were also excluded from fitting, because of large response uncertainties 
and poor statistics, respectively. Because Cyg X-1 was bright, we added a systematic errors of 1%, 
representing the NXB reproducibility (Fukazawa et al. 2009), to all bins. 

In order to grasp a spectral shape and reconfirm the analysis procedure in Paper I under a con- 
dition of using only the HXD data, we fitted the time-averaged HXD spectrum of the first observation 
with a single power-law model. The interstellar absorption was expressed by wabs in xspec, and its 
column density was fixed at A^h = 6 x 10^^ cm~^ after Paper I since its effects are minor in the HXD 
energy band. The best fit photon index was obtained as F = 1 .78. However, the fit was far from being 
acceptable, with x^ /^= 4262/132, because the data are as shown in figure 5 (b) more convex than 
the model. This agrees with a number of past reports, including Paper I. Therefore, we substituted a 
cutoff power-law model for the single power-law model. The fit was improved to x^/'^=170/131, with 
F =1.39 and the exponential cutoff energy _E'cut=162 keV. However, the fit is not yet fully acceptable, 
because of some excess in the 20-40 keV range as seen in figure 5 (c). 

To account for this structure, we introduced a reflection component after Paper I. The utilized 
model is again pexrav (Magdziarz & Zdziarski 1995), which takes into account reflection from cold 
matter and a spectral cutoff. Like in Paper I, we fixed the inclination angle of the reflection plane at 
45°. As shown in figure 5 (d), the fit became fully acceptable with F =1.44, £'cut=158 keV, a reflection 
solid angle l]/27r=0.14 and xVi^=148/130. 

Now that an acceptable fit was achieved in the HXD range by an empirical cutoff power- 



law model including reflection, the next step is to employ the more meaningful description of the 
spectrum, namely unsaturated inverse Comptonization processes (Paper I). For this purpose, we uti- 
lized compPS model (Poutanen & Svensson 1996), which takes into account not only relativistic 
Comptonization correctly but also reflection by a cold matter. Free parameters of the model in our 
analysis are electron temperature kT^, optical depth r, reflection solid angle Vt/2ii and normalization. 
The seed photon temperature was assumed to be 0.2 keV (Paper I), since this parameter cannot be 
determined by the present data which are limited to above 10 keV. We assumed a spherical geometry 
in the spectral analysis. Then, as listed in table 2, the CompPS fit to Observation 1 has become fully 
acceptable, but has yielded the parameters which are not fully consistent with those of the higher r 
component derived in Paper I. The reason of these deviations is discussed in section 5. 

We next applied this model to the same three observations as in figure 2, and obtained the 
results shown in figure 6 and table 2. Since the three fits were all successful, we further applied the 
same model to the remaining observations, and again obtained acceptable fits in all of them (table 2). 
In figure 7, the obtained parameters are expressed as a function of Casm- From figure 7 (a), the HXD 
flux is seen to increase with Casm, at least up to Casm ~ 20 cts s~^ The values of kT^ = 70-100 keV 
and T = 1-1.5 are consistent with those reported in previous works (e.g. Shapiro et al. 1976; Sunyaev 

6 Trumper 1979; Sunyaev & Titarchuk 1980; Gierlihski et al. 1997; Di Salvo et al. 2001; Frontera et 
al. 2001a; Zdziarski & Gierlihski 2004; Ibragimov et al. 2005). 

In figure 7, neither kT^ nor r exhibits significant dependence on Casm- However, as expected 
from the spectral slope change observed in figure 6, a clear anti-correlation is revealed in figure 

7 (b) between Casm and the y-parameter, y = ATkTc/{mt.(?), where k is the Boltzmann constant, 
rrte is the electron mass, and c is the light velocity. Moderate values of reflection, VL/2n: = 0.20- 
0.35, are consistent with previous works, considering that Suzaku covered near the hardest end of 
the LHS (figure 1). As already suggested by figure 6, n/2n depends positively on Casm- Similar 
positive relations between T and Q/2n have been observed not only from Cyg X-1 (Gilfanov et al. 
1999; Ibragimov et al. 2005; Titarchuk et al. 2007), but also from other BHBs and Seyfert galaxies 
(Zdziarski et al. 1999). 

4. Timing Analysis 

4.1. Power Spectral Density 

By using powspec in XRONOS, we computed power spectral densities (PSDs) of all observa- 
tions, separately on the 10-60 keV PIN events and on the 60-200 keV GSO events. The NXB events, 
which account for ~ 50% of the GSO data and -- 2% of the PIN data, were subtracted. PSDs of Cyg 
X-1 in the LHS are known to have a nearly constant power density in frequencies below z/b ~ 0. 1 Hz, 
and becomes "red" above i^b (Miyamoto et al. 1989; Negoro et al. 2001; Titarchuk et al. 2006). This 
i^b is referred to as a "break frequency" and thought to have some relations to variation time scales of 
BHBs, perhaps a viscous time scale (Done et al, 2007; Ingram & Done, 2011). Therefore, our PSDs 



Table 2. Pai'ameters obtained by fitting tlie HXD spectra witli a single CompPS model . 

No. 10-400 keV flux y kT^ r Q./2'k xlii^) & d.o.f. Casm 

(10"* erg cm"^ s"^) (keV) (cts s~^) 

~ 4.74l«;- O.mt^Z 80.1+rj 1.39+g:g? 0.20«:»^3 

2 4.48l°:°« 0.9361°:°°^ 76.41^;? 1.57t°:S? 0.18«:^? 

o Q 9Q+0.03 nQS'J+O-OOS 8^9+3-7 i Al+O-O'T niQ+0-03 

J 0.zy_y Q4 U.;?JJ_Q yog OJ.Z._34 i-^'-0.07 "•^"-0.02 

4 2.09lH^ 0.898t°:«l* 108.6t^«3^ i-Oet^il 0.18«:«^3 

c T 77+0.05 n SQQ+0007 8^ /t+4-8 loo+O.OS r> 7^+0.02 

-J ■-'■''_o.04 "-'■°^"-0.009 °"-^-3.4 ^-^^-0.08 *^-^^-0.02 

6 4.27in? 0.918+°:««^ 80.ltl^ 1.47«:S^ 0.18t°:«^ 

7 4 10+007 r)qf,^+0.012 QS fi+6-5 199+0.10 fi 9^7+0.02 

8 3.681HI 0.917t°:«?« 85.61-^ 1.37t°:i« 0.21«:«^3 

9 3.861°:°! 0.906«:°?« 88.7+-^ 1.31«:?^ 0.22+°:°^ 

10 ,7, +0.06 0884+"°^* QQ7+6-7 1 IC5+O.12 f,99+0.02 
lU J-'i-O.OS "•°°^-0.015 ='''•'-7.9 i-i-^-O.ll "-^^-0.03 

11 4 11+0.05 864+°°i2 92 5+^-5 121+°°'' 27+°°2 

^•^ ^-^^-0.05 "-'■"'^^-O.Oll ^^■•'-6.5 i-^^-O.ll "■^'-0.03 

12 3.49l°:°I 0.820t°:°lJ 93.ll^:« 1.20t°:f8 0.24t°:°^3 

13 4.011°:°! o.8i3t°:°?« 75.2I*:? i.38t°:g^ o.25t°:°^ 

14 3.32l°:°! 0.755«:r8 69.61^^ 1.40tg:°^ 0.27t°:°^3 

15 4.00lg:°^ 0.684«:°1^ 83.7l«:^ 1.05t°:l° 0.32t°:°^3 

16 3 51+0.06 0759+°°" 80 1+1°* ,19+0.14 097+0.04 

10 J.Jl_QQg '^•'-'^-0.023 °"-i-6.2 ^-^^-0.15 "-^'-0.04 

17 3.151°:°^ 0.8621°:°}° 85.21^:1 1.301°:}! 0.201°:°^ 

18 3.421°:°^ 0.8861°:°}? 89.0l«:g 1.28l°:°« 0.211°:°^ 

19 2.751°:°^ 0.8821°:°}° loi.ol}!}:^ 1.091°:?^ o.iit^ol 

20 2.72l°:°! 0.869lg:°^« 109.2l^^,« 1.02l°:}^ 0.26l°:°^3 

21 4.52l°:°I 0.8971°:°°? ^l-^tH 1.42l°:°« 0.271°:°^ 

22 3.74l°:°«3 0.8681°:°}^ 86.511? 1-291°:}° 0.20l°:« 

23 3.821H? 0.897l°:r, 82.7ll« 1.38l°:}J 0.2ll°:°^ 

24 3.80lH^ 0.9051°:°?^ 84.711^, 1.36l°:°^9 0.20l°:°J 

25 5.45lg:}i 0.847lg:g}? 86.91^3 1.26l°o:}^ 0.301^^ 

cover from lO^'^ to 10 Hz, which fully contains the expected Vh- Considering these, we employed 
a binning time of 0.05 s and a calculation interval of 8192 (2^^) bins. Considering our exposure ef- 
ficiency of ~ 50%, any time interval which has less than 4096 bins was discarded. The PSDs are 
normalized so as to directly reflect the fractional rms variation, and are displayed after subtracting the 
expected white noise arising from the counting statistics. 

Figure 8 shows PSDs multiplied by the frequency of the three representative observations 
which were selected in subsection 2.3. All of them exhibit shapes typical of this BHB, while the 
value of z/b and the power density below z/b are both seen to depend on the observation. Interestingly, 
the PSD behavior above z/b is relatively the same, both in the slope and intensity, among the three data 
sets. 

To quantify these PSDs, let us fit them with an empirical formula as 

8 



1.07(130) 


28.6±0.3 


1.00(134) 


21.0±0.6 


1.15(134) 


14.9±0.4 


0.74(134) 


13.5±1.0 


1.09(134) 


15.3±0.2 


1.08(134) 


21.7±0.3 


1.17(134) 


21.8±0.4 


1.20(134) 


17.6±0.3 


1.11(134) 


19.7±0.4 


1.07(134) 


20.8±0.4 


1.06(134) 


22.4±0.4 


1.00(134) 


23.6±0.3 


0.91(134) 


31.8±0.4 


1.38(134) 


33.7±0.4 


1.16(134) 


45.2±0.4 


1.04(125) 


41.8±0.5 


1.21(133) 


20.0±0.3 


1.12(134) 


16.1±0.5 


0.94(129) 


16.6±0.4 


1.26(133) 


15.1±0.6 


1.11(134) 


20.2±0.8 


0.97(132) 


21.8±0.4 


0.94(134) 


25.1±0.4 


1.18(134) 


22.8±0.4 


1.00(118) 


40.1±1.2 



where v is the frequency, while the free parameters, a and c, are normalization and the power-law 
index after the break, respectively. According to RXTE PCA results, a PSD of Cyg X-1 in the LHS 
has roughly two break frequencies (Nowak et al. 1999; Pottschmidt et al. 2003). Our v\, in equation 
1, at 0.03-1 Hz in figure 8, corresponds to their lower break frequency, while the effects of the higher 
break frequency are approximately accounted for by the parameter c. Only the PIN PSDs were fitted, 
because the GSO data give much poorly determined PSDs. 

We obtained acceptable fits with equation 1 to all the 25 HXD-PIN PSDs, when we assume 
systematic errors of at least 3%. Figure 9 summarizes, as a function of Casm^ the obtained PSD 
parameters, where the low frequency power in panel (b) was calculated by integrating P{y) from 
10~^ Hz to 10^^ Hz. As Casm tripled, the value of v^, increased from 0.03 to 0.3 Hz, and the variation 
power over 10~^-10^^ Hz decreased by a factor of ~ 30 (or the fractional variation decreased by a 
factor of ~ 5). These mean that the time scale of variability shortens, and the variation amplitude in 
low frequencies decreases as Casm increases (or the energy spectrum softens). 

Since figure 9(a) was integrated only over the very low frequency range (to avoid the changes 
in z/b), the implied variation amplitude is relatively low, i.e., 2-8% of the mean. In order to more 
properly estimate the variation amplitude and its energy dependence, we divided the three repre- 
sentative HXD data into 6 energy bands, and calculated fractional rms variations by integrating the 
corresponding PSDs over the entire frequency range of lO^'^-lO Hz. The Poisson noise contribution 
was subtracted. As figure 10 shows, the obtained fractional rms variations are 20-40%, with a nega- 
tive dependence on Casm, but little (or slightly negative) dependence on the energy. These properties 
are consistent with previous measurements (e.g., Gierlihski et al. 2010) 

4.2. Auto Correlation Function 

We computed auto correlation functions (ACFs) of the 10-60 keV PIN data and the 60-200 
keV GSO data, as our second step for characterization of the time variability. Although an ACF 
and the corresponding PSD, being Fourier conjugate to each other, carry essentially the identical 
information, the former works in time domain whereas the latter in frequency domain. Therefore, we 
may benefit by studying both of them, and examining them for their mutual consistency. The utilized 
tool is autocor in XRONOS, with a binning time 0.1 s and 4096 bins in each computation interval. 
In other words, we computed the ACFs every 409.6 s with a time resolution of 0.1 s, and took their 
average. We chose to normalize the ACF values by the variances and the number of bins per interval, 
and subtracted the Poisson noise contribution. 

Figure 1 1 displays the ACFs calculated with the PIN and the GSO data for the three repre- 
sentative observations. Thus, the ACFs monotonically decrease on typical time-lag scales of 0.2-1 s, 
which become shorter as Casm get higher. This behavior is in qualitative agreement with that of v\,, 
and the time-lag scale is consistent with the corresponding (27ri^b)~^- 
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We tried quantifying these ACFs by calculating their characteristic time scales using their 
standard deviation from the zero lag, namely, 



a 



Y..;ACF{U)*t^ 



(2) 



\ E. AC Fit,) 

where ti is the time lag at the i-th bin. Hereafter, we denote a for the PIN ACF as crpiN, and that for 
the GSO ACF as ctgso- For the PIN data, the summation was taken from zero lag to the point where 
correlation value first reaches 0.01. In calculating (TgsOj the summation was taken over the same time 
lag range as the corresponding PIN data. Figure 12 (a) shows the obtained crpiN against i^b- The tight 
anti-correlation between the two quantities confirms the expected consistency between the PSDs and 
the ACFs. 

Since an ACF is generally less subject to the Poisson noise than the corresponding PSD, the 
GSO ACFs, such as shown in figure 1 1 (green), allow us to study the reported energy dependence of 
the Cyg X-1 variability (Miyamoto et al. 1989; Maccarone et al. 2000). Actually, figure 11 suggests 
that ctgso is somewhat smaller than ctpin- In order to examine this possibility, o-piN is plotted against 
(Tgso in figure 12 (b). As suggested by the three representatives, (Tgso is thus confirmed to be always 
smaller than o"pin by ~ 20 % in our data sets. 

4.3. Cross Correlation Function 

Previous GINGA observations of Cyg X-1 discovered interesting "phase-lag" phenomena 
among variations in different energy bands (Miyamoto et al. 1989). As a brief examination of such 
effects in the present data which span a much broader energy band than the Ginga data (typically 
2-30 keV), we computed, using crosscor in XRONOS, cross correlation functions (CCFs) between 
the 10-60 keV PIN data and the 60-200 keV GSO data, with the same conditions as employed in 
the ACF analysis. More detailed cross spectrum analysis is a subject of our future work. Each CCF 
was normalized to the geometric mean of the underlying two ACFs, from which the Poisson noise 
contribution was subtracted as in subsection 4.2. 

The obtained CCFs of the three representatives are shown in figure 13. Unlike the ACFs, their 
peak values become less than 1 if variations in the two energy bands are not perfectly correlated with 
each other. Let us examine these CCFs for their lower-to-higher-order moments. As to the first- 
order property, the three CCFs in figure 13 are all peaked at zero time lag within the employed time 
resolution of 0.1 s. This is consistent with previous works in energies below ~ 100 keV, including an 
RXTE upper limit of 2 ms on time lags between 2-5 and 24-40 keV signals (Maccarone et al. 2000). 

Next, from the second-order view, we confirm in figure 1 1 that the characteristic time scale 
of CCFs gets significantly shorter toward higher values of Casm, as readily expected from the ACF 
analysis. What becomes possible with CCFs is asymmetry against the time reversal, namely the third- 
order moment, which is related to phase-lag properties between the two energy bands. In fact, we find 
in figure 13 that the characteristic time scale of correlation toward the positive time lag is longer than 
toward the negative time lag. This means that higher energy photons have some delayed components 
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against softer ones, in qualitative argument with the reported hard-phase-lags (Miyamoto et al. 1989; 
Nowak et al. 1999; Maccarone et al. 2000; Kotov et al. 2001; Pottschmidt et al 2003; Arevalo & 
Uttley 2005). 

To study this property more quantitatively, we calculated skewness as, 
^_ E.CCFiU)*ii-ur 

E.CCFiU) * lE^CCFiU) * (t-ml 
where i is the centroid of the CCF distribution. In every observation, i was calculated as ~ 0.1 s 
or smaller. The range employed for this summation is the standard deviation from the center of each 
CCF since the calculated skewness, the third-order quantity of time, could have too large uncertainties 
if we continued the calculation until zero correlation. By the definition of our CCF time lag, S > 
means hard-phase lag. The results of these calculations, given in figure 14, reveal that the skewness 
of every CCF is larger than 0, so there are hard-lags in all observations. 

5. Discussion 

We have analyzed the Suzaku HXD data of Cyg X-1, acquired on 25 occasions from 2005 
September through 2009 December, and studied the 10-400 keV spectral properties as well as the 
10~'^-10 Hz timing behavior. The source resided in the LHS throughout, accompanied by factor ~ 3 
changes in the 1.5-12 keV RXTE ASM intensity, Casm, which was used as an indicator of the soft 
X-ray flux and of the mass accretion rate. Among these data sets, the first one was already analyzed 
in Paper I. Thanks to the high performance of the HXD including its GSO scintillators, the present 
study has provided considerably improved hard X-ray information on this leading BHB, particularly 
in energies above ~ 100 keV, compared with previous studies (e.g. Miyamoto et al. 1989; Gierlihski 
et al. 1997; Gilfanov et al. 1999; Nowak et al. 1999; Maccarone et al. 2000; Di Salvo et al. 2001; 
Frontera et al. 2001a; Pottschmidt et al. 2003; Zdziarski & Gierlihski 2004; Ibragimov et al. 2005; 
Titarchuk et al. 2007; Gierlihski et al. 2010). 

5.1. Results and interpretation from spectral analyses 

The 25 HXD spectra covering 10-400 keV were all reproduced successfully by a single 
compPS model, incorporating reflection from a cold matter. Therefore, as envisaged repeatedly 
in the literature (section 1) including Paper I in particular, the hard X-ray production in Cyg X-1 
can be interpreted in the disk-corona framework; a disk provides soft X-ray seed photons, which are 
Comptonized by a hot corona, and the produced hard X-ray photons are partially reflected by the disk. 

In terms of these compPS fits, the 10-400 keV energy flux in the 25 observations was calcu- 
lated as (2-5) X 10"^ erg cm~^ s~^ [figure 7 (a)]. Through very crude analysis of some representative 
of the XIS data, we estimated the 0. 1-10 keV flux in our observations as (2-4) x 10~^ erg cm~^ s~^, 
after removing photoelectric absorption. Therefore, the unabsorbed flux in the 0.1-400 keV range, 
where most of the emission is contained (cf. figure 10 in Paper I), is estimated as (4-8) xlO~^ erg 
cm~^ s~^ in our observations. Assuming isotropic emission and a distance of 2 kpc (Ziolkowski 
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2005), these values correspond to 0.8-1.6% of the Eddington luminosity for a black hole of 16 solar 
masses (Caballero-Nieves et al. 2009). It would be interesting to compare the state of Cyg X-1 be- 
tween our data sets and a large amount of pointed observations conducted by RXTE (Gierlihski et al. 
2010). In terms of the "state parameter" defined in Gierlihski et al. (2010), the 25 Suzaku data are 
found in the range of 0-0.25, which is very similar to the LHS covered by the RXTE observations. 
The present results thus provide unified broad-band information on this leading BHB over a typical, 
and reasonably wide, flux range in the LHS. 

The wide-band spectroscopy with Suzaku has successfully decoupled the intrinsic thermal 
cutoff in the continuum and the reflection hump, and gave the reflection solid angle r2/27r as 0.20- 
0.35. This parameters was previously obtained by Gilfanov et al. (1999) as 0.4-0.6 with RXTE, by 
Ibragimov et al. (2005) as 0.2-0.5 with Ginga, OSSE and RXTE, while as 0.1-0.2 by Gierlihski et al. 
(2010) with RXTE. Thus, the absolute values can vary depending on the model and the band width 
utilized for analysis. The cool disk is therefore thought to penetrate halfway into the corona (Paper 
I), because these values of Vl/2t: would be too small if the disk reached the innermost stable circular 
orbit, while too large if the disk were detached from the corona. 

The compPS fits have revealed some interesting dependences of the Compton parameters on 
Casm (figure 7). Although neither Tp nor r exhibits noticeable systematic behaviors, their product, y, 
clearly decreases as Casm increases [figure 7(b)]. Therefore, the effects of Comptonization on indi- 
vidual photons are considered to diminish as the soft X-ray flux becomes higher. Then, the positive 
dependence of the HXD flux on Casm [figure 7(a)] must be a result of higher seed photon inputs 
to the Compton corona towards higher mass accretion rates. Further considering the positive r2/27r 
vs Casm correlations as noted above,the disk-corona geometry must evolve towards higher accretion 
rates in such a way that the energy transfer from the corona to the disk photons decreases, while the 
seed photon input and the disk reflection both increase. After discussing the time variations, we come 
back to this issue in subsection 5.3. 

Our single-zone modeling is admittedly rather approximate, because the corona, in reality, is 
likely to be inhomogeneous; a full-band spectral modeling including the XIS data would require mul- 
tiple Compton components with different ^/-parameters (Paper I). Actually, our results on Observation 
1 are somewhat discrepant from those of the harder Compton component which comprises the double 
Comptonization model employed in Paper I. For example, our y-parameter, 0.87 (table 2), is smaller 
(softer) than that of the harder Compton component in Paper I, y = l.lStaoa, because in that model- 
ing some low-energy signals were taken up by the softer Compton component. The reflection solid 
angle we obtained (fi/27r = 0.20; table 2) could also be smaller than that of Paper I {0.1271 = 0.4+o 3), 
probably because a slightly concave double Compton continuum requires a stronger reflection to 
compensate for it. Nevertheless, we consider our single-zone approach to be meaningful, because the 
single compPS parameters derived here from the HXD data are regarded as some kind of averages 
over the suggested coronal inhomogeneities. 

Incidentally, we obtained somewhat lower values of Tc than in Paper I. This effect is mostly 
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due to the updated GSO data reprocessing (Yamada et al. 201 1, in press), rather than to the modeling 
difference. 

5.2. Results and interpretation from timing analysis 

The HXD onboard Suzaku, with a relatively small effective area (e.g., ~ 270 cm^ at 100 keV; 
figure 10 of Takahashi et al. 2007), is not necessarily optimized for studies of fast aperiodic variations. 
Nevertheless, its extremely low background and wide energy band have allowed us to obtain some 
improved information on the fast flickering of Cyg X-1, particularly in energies above ~ 100 keV 
where previous hard X-ray instruments were limited by background. 

Our timing analysis gave 10~^-10 Hz PSDs (figure 8) which are typical of this object; a 
relatively flat power distribution and a red-noise behavior, below and above a rather well defined break 
frequency i^b^ respectively. The values of u\, we found, 0.03-0.3 Hz [figure 9 (b)], are consistent with 
the range of 0.03-0.5 Hz which Pottschmidt et al. (2003) derived from the 2-13 keV RXTE PCA 
data wherein Casm varied over 15-80 cts s^^. Our ACF studies have been consistent with the PSD 
analyses [figure 11, figure 12 (a)]. 

Figure 10 shows that the total fractional rms variations, integrated over 10~^-10 Hz, were 
found in the range of 20-40%, which is typical of this object. In addition, the fractional rms variation 
was found to be approximately energy independent up to 100-200 keV. This reconfirms our basic 
interpretation (section 1) that the X-ray flickering reflects variations in the seed photon input to the 
corona. 

As already mentioned in subsection 4.1, one important result from our timing studies is the 
prominent evolution of the PSDs (figure 8, figure 9): towards higher values of Casm^ ^h increases and 
the low-frequency power decreases, while the PSD above v\, is kept rather unchanged. The change 
in z/b is also supported by the ACF behavior in figure 1 1 . Then, without invoking particular models 
of variability, we may argue in the following way. The observed aperiodic flux changes may be 
regarded as a superposition of some local (or propagating) fluctuations, occurring in the corona (or in 
the coronal covering fraction over the cool disk) which is thought to span a broad radius range of a 
few to ~ 100 gravitational radii (Paper I). Furthermore, it is plausible to think that slower (or faster) 
variations are produced mainly in outer (or inner) regions of the corona, where various physical 
phenomena generally have longer (or shorter) time scales. Then, one of the most natural ideas would 
be to presume that an increase in the mass accretion rate causes the outer coronal boundary to shrink 
inwards, and hence the slowest variations produced in the outermost region disappear. This can 
consistently explain the observed three effects; the increase in z/b> the decrease in low-frequency 
power, and the unchanged high-frequency power. 

To be somewhat more specific, let us assume, for example, that the variation time scale is re- 
lated to viscous time scales in an standard "alpha" disk (Shakura & Sunyaev 1973), which is described 
as 

U,(xa-\r/Hf^K-' . (4) 
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Here, a. is the viscosity parameter, r is the radius, H is the scale height of the accretion disk or the 
corona, and fix is the Keplerian angular velocity. Supposing that a and rjE are both kept unchanged, 
then coronal size changes by a factor of ~ 5 are sufficient to explain the observed change in i/b by an 
order of magnitude among the present data sets. 

53. The inferred accretion geometry 

We may now combine the spectral view developed in subsection 5.1 with the variability con- 
sideration in subsection 5.2. In the former subsection, we observed the three spectral evolutionary 
changes as Casm increased; the increase in the seed photon input, the decrease in y, and the increase 
in the reflection solid angle. These results were interpreted in terms of a view developed in Paper I, 
namely, partial overlap between the disk and the corona. In the latter subsection, an increase in the 
mass accretion rate was found to cause three effects, namely, the increase in v\,, the increase in the 
low-frequnecy power and the rather unchanged high-frequency power density. We interpreted these 
results that the outer boundary of the corona, to be denoted i?c, shrank (e.g., by a factor of 5 from the 
v\, change) as the mass accretion rate increased. 

When i?c decreases towards higher accretion rates, the coronal height H is expected to also 
reduce in rough proportion to R^., because the corona is considered to be approximately spherical. 
As y is approximately proportional to ricTcH (with n^ the average electron density), the observed 
concurrent decrease in y can be explained if the coronal electron pressure riekTc somehow remains 
approximately constant, or if it even decreases, e.g., under increased Compton cooling. 

Let us consider the behavior of the innermost disk radius i?in, in order to explain the remaining 
two spectral evolutions. One of them is a need for a larger number of seed photons to be supplied to 
the corona at higher accretion rates. This clearly implies a larger overlap between the corona and 
disk, thus requiring Ri^ to decrease by a larger factor than Re. This view is based on a conclusion 
of Paper I that the cool disk is truncated at ~ 15 Rq. Then, the remaining property, namely the 
increased reflection, follows as a natural consequence, because a larger fraction of the corona becomes 
penetrated by the cool disk. The invoked larger disk-corona overlap will also decrease T^ of the corona 
through increased Compton cooling, and will contribute to the required reduction in y. Thus, the 
observed spectral evolution, as well as those in the variation time scales, can be explained consistently 
and very naturally by assuming that the disk intrudes more deeply into the corona as Re decreases. 

5.4. Implications of the energy dependence of time variability 

Finally, let us discuss implications of the subtle but intriguing two energy-dependent effects 
in fast variations, as revealed or reconfirmed utilizing the genuine broad-band capability of the HXD. 
One of them is the relation of ctgso < o"pin> found in all observations by the ACF analysis (figure 12). 
In figure 8, this effect may be visible as a slightly flatter slope above ~ 0.3 Hz in the GSO PSD than in 
the PIN PSD. In other words, harder photons vary systematically on shorter time scales as observed 
previously (e.g., Maccarone et al. 2000; Pottschmidt et al. 2003), and the effect persists at different 
mass accretion rates. 
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In the present disk-corona picture, the faster variability in harder photons can be explained 
most naturally, if the multi-?/ property revealed in Paper I takes a form of negative radial gradient in y. 
Then, harder photons are expected to vary on faster time scales, because they will be produced more 
efficiently in regions with larger y, namely inner regions, where the variability time scales should be 
shorter. Though still speculative at present, a particularly attractive possibility is that the outer coronal 
region overlapping with the disk has lower y value(s) under photon cooling, while the non-overlapping 
inner region has higher y. 

The other energy-dependent effect is the non-zero skewness seen in figure 14, which is essen- 
tially equivalent to the CCF asymmetry: while variations in the 60-200 keV range exhibit no direct 
time lags within 0.1 s from those in 10-60 keV, their cross correlation is stronger when the harder 
signals are delayed (figure 13). These results extend, up to ~ 200 keV, the previously reported "hard- 
phase lags", namely, the effect that a cross spectrum between softer and harder signals exhibits phase 
differences that depend much less strongly on the frequency than in the case of a constant time lag 
(e.g., Miyamoto et al. 1989; Nowak et al. 1999; Poutanen 2001; Pottschmidt et al. 2003). Similar 
hard-phase lags have been observed from other BHBs in the LHS (e.g., Miyamoto et al. 1992; van 
der Hooft 1999ab; Nowak et al. 1999). As pointed out in many previous works (e.g., Maccarone et 
al. 2000), this phenomenon in Cyg X-1 and other BHBs is difficult to explain in terms of Compton 
delays, since this would predict (Tgso > ctpin in contradiction to the observations, and would require 
too large a corona to produce phase lags on time scales of ~ 1 s. 

A more reasonable explanation to the hard-phase lags is to assume that fluctuations responsible 
for the X-ray variations propagate inwards through the corona which has a radial gradient in y as 
invoked above (Lyubarskii 1997; Done et al. 2007). Then, harder signals, produced in regions closer 
to the black hole, would be generally delayed from softer ones to be produced in more outer regions. 
Furthermore, the delay time between the two bands would depend on the frequency of variation under 
consideration: the delay time should be shorter for higher-frequency fluctuations, because they must 
originate in regions closer to the black hole, where propagation delays would also be shorter. Then, we 
expect the energy dependence to be of phase-lag type than a constant-time-lag type. These ideas were 
already presented by Nowak et al. (1999), although their main focus was on the Compton scattering. 

In this way, the two energy dependent effects of time variability have reinforced the mxxlii-y 
picture derived in Paper I. While Paper I assumed a condition wherein domains with high and low 
values of r coexist at each radius, it has become more likely that a large-scale gradient in y is also 
present. 

6. Conclusion 

Analyzing the Suzaku HXD data of Cyg X-1 acquired on 25 observations over 2005 October 
through 2009 December, we have obtained the following results. 

1. Cyg X-1 was in the LHS on these occasions, accompanied by a factor 3 changes in Casm- The 
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estimated 0.1-400 keV unabsorbed flux varied over (4-8) x 10~^ erg cm~^ s~^. 

2. The 10-400 keV HXD spectra were all reproduced successfully with a single Comptonization 
model, considering reflection from a cold matter. 

3. As Casm became higher, the HXD flux and the reflection solid angle increased, while the 
Compton ^/-parameter decreased. 

4. When Casm increased by a factor of ~ 3, the PSD break frequency became higher by an order 
of magnitude, and the power integrated over 10^^-10^^ Hz diminished by a factor of ~ 5. 

5. In the ACF analysis, the 60-200 keV signals always showed shorter correlation lengths than the 
10-60 keV ones. 

6. The CCF revealed skewness so that the 60-200 keV signals were delayed from those in the 10-60 
keV, and the effects became more prominent towards higher mass accretion rate. 

In order to consistently explain items 3 and 4 above, we have developed an interpretation that a 
gradual increase in the mass accretion rate causes the corona to diminish in size, while the cool disk to 
intrude more deeply into the corona. Property 5 and 6 listed above may be explained by considering a 
radial gradient in the corona, so that harder photons are produced in closer vicinity of the black hole. 
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Fig. 1. A scatter plot between the 1.5-12.0 keV count rate and the 5.0-12.0 keV vs 1.5-3.0 keV hardness ratio of 
Cyg X-1, obtained with the RXTE ASM over a period of 1996 December to 2010 September Each gray dot is an 
integration over one day. Black points indicate epochs of the 25 Suzaku observations. 
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Fig. 2. Background-subtracted 64-s bin light curves of Cyg X-1 with HXD-PIN (10-60 keV;top panels) and 
HXD-GSO (60-200 keV;bottom panels) from three representative observations. Left two panels are for Observation 
3, middle two panels for Observation 10, while the right two panels Observation 15. 
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Fig, 3. Time-averaged and response-inclusive spectra of HXD-PIN (black) and HXD-GSO (red) from the three 
representative observations, show after subtracting the simulated background data denoted as PIN_bkgd (blue) and 
GSO_bkgd (magenta), respectively. Also shown in magenta are 3% and 5% levels of the GSO background. Panels 
(a), (b) and (c) correspond to the three epochs selected in figure 2. 
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Fig. 4. (a) The count rate of PIN and GSO, and (b) the HXD hardness ratio (see text) of Cyg X-1 plotted against 

Casm- 



10 
1 
M 0.1 
I 0.01 
0.001 

1 

0.7 
.2 0.4 
I ^-^ 

T3 1 

O i 

I 0.7 

Q 1-3 

1 

0.7 



Observation 1 



(a) pexrav 



^.^^mfii ----. 



(b) power-law (r=1.82) 



(c) cutoffpl 






^; 



..j...AI^L,lu,Llli 



(d) pexrav 



A 



^ 



A 



10 20 50 100 200 400 

Energy (keV) 

Fig. 5. (a) The response-folded PIN and GSO spectra from the first observation, compared with the best-fit pexrav 
model. Panels (b), (c) and (d) shows the data divided the best-fit powerlaw model, cutoffpl model, and the 
pexrav model (corresponding to panel a), respectively. 
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Fig. 6. The same HXD spectra of Cyg X-1 as in figure 3, presented in the deconvolved lyFv form. The best-fit model 
and the fit residuals are also shown, where the Compton component and the reflection component are displayed 
separately. 
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Fig. 7. The Compton (CompPS) parameters from the 25 observations of Cyg X-1, plotted against Casm- (a) The 
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Fig. 8. Power spectral densities (PSDs), multiplied by the frequency ly, of the PIN data (panel a) and the GSO data 
(panel b) from the three representative observations. Red shows Observation 3, green Observation 10 and blue 
Observation 15. 
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Fig. 9. Correlation between Casm and (a) power integrated over 10 •^-10 ^ Hz, and (b) the break frequency. 
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Fig. 10. The fractional rms variation spectra of the three representative observations. 
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Fig. 11. The PIN (red) and the GSO (green) ACFs of the three representative observations, calculated in the 10-60 
keV and 60-200 keV range, respectively, (a) Observation 3, (b) Observation 10, and (c) Observation 15. 
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Fig. 12. (a) The characteristic time scale of the 10-60 keV PIN ACF ctpin, shown against i^h- (b) Correlation 
between ctpin and ctgso- 
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Fig. 13. CCFs between the 10-60 keV PIN data and the 60-200 keV GSO data from the three representative obser- 
vations (red), with time-inverse ones (green). The points at zero lag are presented with black points. 
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Fig. 14. Dependence of the CCF skewness on Casm- The data point for Observation 16 (Casm = 40.2 cts s ^) is 
omitted, because of too large errors. 
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